Evidence for infrared finite coupling in Sudakov resummation 
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New arguments are presented in favor of the infrared finite coupling approach to power correc- 
tions in the context of Sudakov resummation. The more regular infrared behavior of some peculiar 
combinations of Sudakov anomalous dimensions, free of Landau singularities at large Nf, is pointed 
out. A general conflict between the infrared finite coupling and infrared renormalon approaches to 
power corrections is explained, and a possible resolution is proposed, which makes use of the arbi- 
trariness of the choice of constant terms in the Sudakov exponent. A simple ansatz for a 'universal' 
non-perturbative Sudakov effective coupling at large Nf emerges from these considerations. An 
alternative evidence for an infrared finite perturbative effective coupling in the Drell-Yan process at 
large Nf (albeit at odds with the infrared renormalon argument) is found within the framework of 
Sudakov resummation for eikonal cross sections of Laenen, Sterman and Vogelsang. 
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The notion of an infrared (IR) finite coupling, and the 
related concept of universality, to parametrize power cor- 
rections in QCD has attracted much attention for a long 
time I n the present note I display further evidence 
in favor of this assumption in the more specific frame- 
work of Sudakov resummation, taking the examples of 
scaling violation in deep inelastic scattering (DIS) and in 
the Drell-Yan process (more details can be found in 0). 
From the standard exponentiation formulas in Mellin N- 
space for the DIS structure function (I shall adopt in 
general the notations of 0]), one gets the large N relation 
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is given as a power series m 
= a s (Q 2 )/4ir with A^-mdependent coeffi- 
4C F As(fc 2 ) = A(a s (k 2 )) + dB{a s (k 2 ))/d\nk 2 , 
(the universal 'cusp' anomalous dimension) and 
the standard Sudakov anomalous dimensions 
relevant to DIS, given as power series in a s . Ag(k 2 ) 
shall be referred to as the 'Sudakov effective coupling', 
a physical 'effective charge' j|| and we shall see (in 
the large Nf limit) that it has drastically different 
IR properties then the individual Sudakov anomalous 
dimensions it is composed of. 

The Borel transform B\As](u) of the Sudakov effective 
coupling has been computed |f| at large Nf . It is defined 
by A s {k 2 ) = j- Q r o °° duexp(-u\nk 2 /A 2 )B[A s }(u), where 

(3o = -y C2(G) — ^Nf is the one- loop coefficient of the 
beta function. One gets, with the 'naive non-abelization' 
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recipie: B[A s ](u) = cxp(- du)^^ - ^— , 

where d is a scheme-dependent constant related to the 
renormalization of fcrmion loops: d — —5/3 in the MS 
scheme and d = in the so called 'V-scheme'. In the 
following I shall use the 'V-scheme' for simplicity. I 



note that B[As)(u) is free of renormalon singularities, 
and that for k — A (the Landau pole of the one-loop 
V-scheme coupling) As{k 2 = A 2 ) is finite. This behav- 
ior is in striking contrast with that of the cusp anoma- 
lous dimension 6] which displays wild oscillations, result- 
ing in a completely unphysical behavior around k = A. 
It is actually possible to get an analytic expression for 
As{k 2 ), valid at all k 2 , and one finds that Ag(k 2 ) ap- 
proaches an infinite (and negative) IR fixed point for 
k 2 - 0: A s (k 2 ) ~ + This is not by 

itself an unphysical behavior ( the coupling is causal and 
there is no Landau singularity), except for the negative 
sign in the infrared, which cannot reproduce a vanishing 
Sudakov tail. The other trouble is in the too strongly 
divergent IR behavior, which gives a divergent contri- 
bution to the integral on the right-hand side of eq.Q. 
These facts make plausible the speculation that there ex- 
ists a non-perturbative modification 5A$ (k 2 ) of the cou- 
pling at small scales pj which might turn the infinite IR 
fixed point of perturbative origin into a genuinely non- 
perturbative, but softer (eventually finite) fixed point. 

However, there is a potential clash between the IR finite 
coupling and the IR renormalons approaches to power 
corrections (very closely connected to the well-known is- 
sue of 1/Q corrections in Drell-Yan) which can be 
summarized as follows. Consider a typical 'renormalon 

integral' R(Q 2 ,N) = jf^-F(k 2 /Q 2 ,N)A s (k 2 ), 
where F(k 2 /Q 2 ,N) shall be referred to as the 
'Sudakov distribution function', and introduce its 
('RS invariant' 0) Borel representation R(Q 2 ,N) = 
^ J °° duexp(—ulnQ 2 / A 2 )B[R](u, N). The Borel trans- 
form of R is given by B[R](u, N) = B[A s ](u) x F(u,N), 
with F(u, N) = J^F(x 7 N)exp(-uhxx). The factor- 
ized form @|>@ of the expression should be noted. Sup- 
pose now B[Ag](u) has a zero at u = uq, which is not 
shared by B[R](u, N). Then F(u,N) must have a pole 
at u = uq. If uq > 0, this means an IR renormalon, and 
F(x,N) contains an O(x uo ) contribution for x — > 0. If 
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this renormalon is not present in B[R](u, N), the stan- 
dard IR renormalon philosophy would conclude that no 
corresponding power correction is present. On the other 
hand, this power correction is still expected in the IR fi- 
nite coupling approach, where A$(k 2 ) is assumed to have 
a finite IR fixed point and the low energy part of the 
renormalon integral R(Q 2 ,N) is well-defined. Then the 
distribution function F(k 2 /Q 2 ,N) can be expanded in 
powers of k 2 for k 2 — > 0, yielding a non-vanishing power 
correction for each term [!j in its low energy expansion, 
parametrized by a low energy moment of the effective 
coupling As(k 2 ). Even if B[R](u, N) is singular, rather 
than non- vanishing and finite, at u = uq, there is still 
a clash, since the singularity of F(u, A) is necessarily 
stronger than that of B[R] (u, A) in presence of a zero of 
B[Ag](u); thus the IR finite coupling approach will pre- 
dict a more enhanced power correction than indicated by 
the renormalon argument. Two examples of this situa- 
tion at large Nf are: 

i) DIS: there B[A$](u) vanishes for any integer u, u > 3. 

ii) Drell- Yam in this case the Sudakov effective cou- 
pling which occurs in the Q 2 derivative of the Drell- 
Yan cross section is given by 10] 4GfAs,dy (k 2 ) = 
A(a s (k 2 )) + ^dDrjY{a>s(k 2 ))/d\nk 2 , and the correspond- 
ing Borel transform at Nf = oo is 6] (dropping an 
inessential exp(cu) factor which corresponds to a change 

of scale): B[As,dy](u) = rfi+u) r(i/2- M ) > which vanishes 
for positive half-integer u. In both cases, the zeroes are 
absent from the corresponding Sudakov exponent, and 
the IR finite coupling approach will lead to the predic- 
tion of extra power corrections in the exponent (N 3 /Q 6 , 
N 4 /Q 8 , etc.. .in DIS, and N/Q, A 3 /Q 3 , etc.. .in Drell- 
Yan). To make progress, we have to understand better 
the origin of zeroes in B[Ag](u). 

DIS case : putting S(Q 2 ,N) = Q dz zN ~^~ l A S [{1 - 
z)Q 2 ], it is useful to write S(Q 2 , A) as a 'renormalon in- 
tegral', where the Sudakov distribution function is given 
by F(k 2 /Q 2 ,N) = {l-tf/Q 2 ) 1 *- 1 -!. Order by order the 
perturbative series of S(Q 2 , N) contain both O(N ) 'con- 
stant terms' and terms which vanish as A — > oo. A mean- 
ingful simplification is achieved by making use of the fol- 
lowing important scaling property: defining e = Nk 2 /Q 2 , 
and F(k 2 /Q 2 ,N) = G(e, N) = (l-e/A)^- 1 -1, and tak- 
ing the limit A — > oo with e fixed one gets a finite result 
G(e,N) -» G(e,oo) = G(e) with G(e) = exp(-e) - 1. 
Let us now redefine the Sudakov exponent by using 
G(e) as the new distribution function, i.e. S(Q 2 ,N) is 

replaced by S stan (Q 2 ,N) = jf ^-G(Nk 2 /Q 2 )A s (k 2 ) 
where the subscript 'stan' means 'standard'. This step 
is legitimate since, order by order in perturbation the- 
ory, S stan {Q 2 ,N) and S{Q 2 ,N) differ only by terms 
which vanish as N — > oo, and thus share the same In A 
and constant terms. I next assume a similar ansatz for a 
modified Sudakov exponent 
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and show that a unique solution for G new (Nk 2 /Q 2 ) and 
Ag ew (k 2 ) exists, under the condition that S new repro- 
duces all divergent In TV terms, together with an (a priori 
arbitrary) given set of constant terms. The latter are 
usually a subset of the set of all constant terms on the 
right hand side of eq.Q, which are uniquely defined, but 
(arbitrarily) split between those included in the Sudakov 
exponent S(Q 2 , N), and those belonging to H(Q 2 ). This 
statement can be checked order by order in perturbation 
theory. Here I give an all-order proof. I first introduce 
a simplification appropriate to the large N limit. The 
perturbative series of S new (Q 2 ,N) still contains 0{1/N) 
terms at large A. It is useful to discard them, and 
find a Borel representation of the corresponding series 
S^l w (Q 2 , N) which contains only the (same) ln p A and 
O(N ) terms, with all the 0(1/N) terms dropped. One 
can show the looked for Borel transform is given by 



f°° de 

exp(u\nN) / — G new (e) exp(— wine 
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B[S% w }(u,N) = B[Ar w )(u) 
G neu ,(oo) 



(3) 



where G nelu (oo) = —1 is a subtraction term, and is fixed, 
since it determines the leading logs. Eq.© displays the 
general form of the large N Borel transform: 



B[SZJ(u,N)=eMu^N)B(u) 



(4) 



It is clear that the iV-dependent part B{u) of the Borel 
transform is unambiguously determined, and cannot be 
changed without changing the A^-dependent terms, i.e. 
the coefficients of all logs of N, which are fixed. On 
the other hand, the constant terms included into the Su- 
dakov exponent are obtained by setting N = 1 in eq.@ 
or Q, and can be changed only by modifying the func- 
tion C new (u) (with C new (0) = 1). Given B(u) from an 
independent calculation, and given an (a priori arbitrary) 
function C new (u), eq.JJJ and (gj determine in principle 
both BlAcf™]^) = C new (u) and G new (e), which com- 
pletes the proof. 

Large Nf result: Specializing now to large Nf to com- 
pute B{u) (and eventually C new {u)) , the result for the 
Borel transform of dlni^G; 2 , N) / d\xiQ 2 \sDG at finite 
N can be given in the 'massive', 'single dressed gluon' 
(SDG) dispersive Minkowskian formalism ,pj as 



B[d In F 2 (Q 2 ,N)/d\n Q 2 ]sdg {u, N) = AC F f - 
f°° dx ■■ 

x / — Fsdg (x, N) exp(-u In x) , 
Jo x 



(5) 



^G new (Nk 2 /Q 2 )A^ ew (k 2 ), (2) 



where the 'characteristic function' !Fsdg{x, N) has been 
computed in Q (x = X 2 /Q 2 where A is the 'gluon mass'). 
I have checked that here too a similar scaling prop- 
erty holds at large N, namely, putting Gsdg(Ui A) = 
T SDG (x, A) with y = Nx = NX 2 /Q 2 , one gets for A -> 
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oo at fixed y: Qsdg{v,N) -> GsDG(y,°°) = Gsdg(v), 
with 



Qsdg{v) = -1 + cxp(-y) - -y exp(-y) 

~y r(0,y) + ± y 2 T(0,y), (6) 

where T(0, y) is the incomplete gamma func- 
tion and implies J °° ^-Qsdg{u) exp(-Mlny) = 



for AT 



i?[dlnF 2 (Q^,7V)/dlnQ 2 ]^ DG ( U ,7V)=4Cf 

ITU 

exp(ulniV) / — Q S dg{v) exp(-u In y) + rs - PG ^'"- ) 
Jo V u 



where Tsdg^) = (^t) 2 jr^KT^m is again a 
subtraction function which takes into account all 
constant terms of the SDG result. The latter are 
indeed obtained by setting N = 1 in eq.J7|). Thus, 
if one wants to select as input an arbitrary subset of 
constant terms to be included into a new asymptotic 
Sudakov exponent S^l w (Q 2 , N), one should simply 
change the subtraction function Tsdg(u), namely 
define a new input B[d\nF 2 (Q 2 ,N)/d\nQ 2 } a s s ] ^ w 
by eq.0, with T SDG (u) -> T new (u) where T new (u) 
(with T new (0) = 1) takes into account the new 
set of constant terms. Having no 0(1/ N) terms, 
one can now identify 4:C F B[S^ s ew ](u, N) with 
B[d\nF 2 (Q 2 ,N)/d\nQ 2 } a s % n £ w (u,N), which yields 
the two master equations 

B[A2 ew ](u) / -G new (e)exp(-u\ne) = B{u) 
Jo e 

sin7ru f°° dy - . 

— QsDG{y)exp(-uhiy), (8) 



and 
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B[A n s ™](u) = C new {u) 



sin 7TU 



,{u). (9) 



gives i?[A^ e ]( W ) 



/?0 



, hence 

— arctan(r/7rj 

2 7T 



with i = lnfc 2 /A 2 . It turns out that A^ mple (k 2 ) is al- 
ready causal and IR finite (with Ag mple (0) = l//3o) 
at the perturbative level, and no non-perturbative mod- 
ification is a priori necessary. However, the first two 
zeroes at u = 1 and u — 2 lead to two leading Zog- 
enhancedO(-Q>i InQ 2 ) and C(-^r hiQ 2 ) power corrections 
at large N in the IR finite coupling framework, as can 
be easily checked since we have in this case G S i mp i e (e) — 
GsD G (e) - i[e(lne + 7f; -3) + e 2 (-ln e - 7B + l) + 0( e 3 )]. 
There is thus, as expected, a discrepancy with the IR 
renormalon expectation, which predicts only 0(^-) and 

O(-^r) power corrections. 

Drcll-Yan case : The analogue of eq.Q for the 
scaling violation of the short distance Drell- 
Yan cross-section is dlnanyiQ 2 , N)/d\nQ 2 = 
4C F [H DY {Q 2 ) + S DY (Q 2 ,N)] + 0(1/N) with 

( 7 )S DY {Q 2 ,N) = J^dz 2^f±A s , DY [(l - z) 2 Q 2 ]. 



Eq. (JHJ) and (J5J) allow to determine both the Sudakov dis- 
tribution function G new (e) and the Borel transform of the 
associated Sudakov effective coupling B[Ag ew ](u), for a 
given input set of constant terms which fix the subtrac- 
tion function T new (u). It is interesting that B[Ag ew ](u) 
is given entirely by the subtraction term. Applying these 
results to the standard exponentiation formula, where 
G(e) is known, and using eq.©, one can in particular 
rederive the previously quoted result for B[As](u). I note 
that the simplest possible ansatz, namely T ' S impie{u) = 1, 



(10) 



Sdy{Q 2 ,N) can again be written as a 'renormalon 
integral', with F(k 2 /Q 2 ,N) -» F DY (k 2 /Q 2 ,N) = 
(1 - k/Q) N - x - 1 and A s {k 2 ) — A s , DY (k 2 ). The 
Sudakov distribution function Fpy (k 2 / Q 2 , N) involves 
both even and odd powers of k at small k. Tak- 
ing the scaling limit A" — > oo with er> Y — Nk/Q 
fixed one thus gets G dy(cdy) = exp(— cdy) — 1- 
Now the work of 12] for eikonal cross sections sug- 
gests to use instead (thus avoiding odd power terms) 

S™™(Q 2 ,N) = jf ^-G^(Nk/Q)A^ Y (k 2 ), with a 
distribution function (I deal with the log-derivative of 
the Drell-Yan cros-section) 



G™ Y (Nk/Q) 



i.e. (j dy (cdy 



d\nQ 2 



[1 



[K (2Nk/Q) + \n{Nk/Q) + lE ] 

(11) 

(x = 2€dy)] , where Kq(x) 



,dKo 
dx 



is the modified Bessel function of the second kind. I note 
that G^fy (ecy) — > —1 for €dy —> oo, consistently with 
the large cdy limit of Gdy(^dy)- At large Nf one gets, 
instead of eq.® (eq.© remains the same) 



f°° de 

B[Af$ Y ]{u) / 2-GSS?(e)exp(-2«ln6) = 
Jo e 

sin7ru f°° dy - . . 

— ysDG,DY{y)exp(-u\ny), (12) 



ttw Jo V 

where y = N 2 X 2 /Q 2 . On the other hand, the large Nf 
calculation of yields the A^-dependent part of the large 
A^ Borel transform 

sin7rii f°° dy « , . 1 T(l — u) 

— Qsdg.dy (y) exp(-ulny) = --^prrr- — r, 

(13) 



tto Jo y 

whereas ea. (|TTJl gives 



/ 2 — G%? (e)exp(-2ulne) = -u[T(-u)} 2 . (14) 
Jo 6 
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From ea. l|12|l one thus derives the large Nf result 

1 sin ttu 



B[A%g Y ]{u) 



r(i + u)r(i - u) nu 



(15) 



hence Ag e $ y 



(k 2 ) = A^ mple (k 2 ). Thus in the framework 
of |l2j | the simplest IR finite perturbative coupling nat- 
urally emerges! The fact that the new Sudakov distri- 
bution function GJfy (e) implies a new Sudakov effec- 
tive coupling is of course one of the main point of the 
present paper. Again, there is a discrepancy [l3| with 
the IR renormalon expectation: although G^y '(cdy) — 
e 2 DY {2\a.eDY + Z^/e — 1) + 0{e\, Y lne^y) involves only 
even powers of tDY for tjyy — ► 0, they are logarithmically 
enhanced. One thus gets Q{^$ InQ 2 ) power corrections 

(p integer) at large N, instead of the 0{^^) corrections 
expected from IR renormalons. 

Large Nf ansatz for a 'universal' Sudakov coupling: the 
question arises whether it is possible to reconcile the IR 
renormalon and IR finite coupling approaches. Actually, 
I should first stress it is not yet clear whether the two ap- 
proaches should be necessarily reconciled. For instance, 
in the DIS case, it could be that the operator product ex- 
pansion (OPE) at large N is consistent with the existence 
of higher order power corrections (N 3 /Q 6 , N 4 /Q s ,...) in 
the exponent. If this turns out to be the case, the IR fi- 
nite coupling approach would be consistent with the OPE 
(and at odds with the IR renormalon prediction) even 
within the standard [j| exponentiation framework. Simi- 
larly, the ansatz Ag mple (k 2 ) might be the correct one in 
the Drell-Yan case although it also contradicts the 
IR renormalon expectation. It is interesting to note that 
eq.@ indicates that zeroes in B[Ag ew ](u) arise from two 
distinct sources: either the 'universal' sin7ra/7ra factor 
{simple zeroes at integer u can come only from there), or 
the 'arbitrary' T new {u) subtraction term (zeroes at half- 
integer u can come only from there). The previous dis- 
cussion suggests that zeroes coming from the 'universal' 
smiru/iTu factor need not be necessarily removed in the 
IR finite coupling approach, at the difference of the more 
'artificial' zeroes (such as those occuring in the standard 
Drell-Yan case) coming from the subtraction term. 
Notwithstanding the above remarks, I shall adopt in the 
following the attitude that the IR renormalon and IR 
finite coupling approaches to power corrections should 
always be made consistent with one another. For this 
purpose, one must remove all zeroes from B[A 7 g ew \(u). 
The mathematically simplest 'universal' ansatz suggested 
by eq.Q is to choose T new (u) = T(l — u), which yields 



7. It is interesting 



B[Ar B ](u) = (^)T(l-u) = r(1+u) . 
to compare this result with the one obtained by elimi- 
nating the half-integer zeroes from the standard Drell- 
Yan result. There the simplest procedure is to define 
B[A% e % Y ](u) = ^p± B[As, DY }(u), which yields again 



the previous ansatz. Given B[Ag ew ](u), eq.JSJ) or fESjl 
determine the Sudakov distribution function G new (e). In 
particular, assuming the above ansatz, one gets in the 
DIS case 



Gnew (^) 



-1 



< e < 1 
£ > 1, 



(16) 



where e = Nk 2 /Q 2 . Thus all power corrections beyond 
the two leading ones are indeed absent from the Sudakov 
exponent, in agreement with the renormalon argument. 
In the Drell-Yan case, one gets instead 



DY { e DY) 



exp(-e|, r ) 



(17) 



i.e. only even powers of eoY — Nk/Q, with no 
logarithmic enhancement for edy — * 0. There is 
compelling numerical evidence that the corresponding 
'universal' Sudakov effective coupling: Ag ew (k 2 ) = 
du exp (— u In fc 2 /A 2 ) l/r(l + u) blows up very 

fast (but remains positive) for k 2 — > 0: Ag ew (k 2 ) ~ 
jr- exp (^jt^J ■ Assuming that this behavior is indeed cor- 
rect, one can again speculate that there exists a non- 
perturbative modification 8A r g ew (k 2 ) of the coupling at 
small scales which will generate a non-perturbative but 
finite IR fixed point (a simple ansatz for SAg ew (k 2 ) has 
been given in |2|). 

To conclude, I have shown that in a number of exam- 
ples at large Nf the Sudakov effective couplings display 
remarkably smooth IR properties, with no Landau sin- 
gularities (a rather unusual occurence in resummed per- 
turbation theory). Generalizing from the present find- 
ings, one can distinguish two main possibilities. Either 
the subtraction term T new (u) has some singularities at 
finite and positive values of u (which must be simple 
poles for the corresponding Sudakov coupling to be free 
of IR renormalons): this is the case in DIS and Drell- 
Yan with the standard choice of the Sudakov distribution 
function; or it has no u > singularities, as the 'simple' 
coupling cci. (|10fl . which emerges naturally in the formal- 
ism of [13] for the Drell-Yan cross-section. In the former 
case, the corresponding Sudakov effective coupling is pre- 
sumably |l4j ] causal, but strongly IR divergent, requir- 
ing a low scale non-perturbative modification to build 
an IR finite coupling. In the latter case, the Sudakov 
coupling is likely (provided T new (u) does not decrease 
faster then an exponential for u — > 00) to be already 
IR finite (and causal) at the perturbative level. Even 
then, there may be extra non-perturbative effects, which 
could also take the form of a non-perturbative modifi- 
cation of the coupling at low scales. I consider these 
facts as hints strongly supporting the IR finite coupling 
hypothesis. The above ansatzes are in conflict with the 
IR renormalon expectations for power corrections. The 
simplest attempt to reconcile the two approaches leads 
to the proposal of a universal perturbative Sudakov cou- 
pling, which again requires a non-perturbative modifi- 
cation to achieve IR finitness. The latter could intro- 
duce some amount of non-universality. The variety of 
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Sudakov distribution functions showing up in the expo- 
nentiation procedure is related to the freedom to select 
an arbitrary set of 'constant terms' in the Sudakov ex- 
ponent. This freedom leaves open the question of the 
determination of the 'correct' Sudakov distribution func- 
tion. There is no obvious way to solve this ambiguity. 
In particular, one can show [2| that the natural option 
of including all O(N ) terms into the Sudakov expo- 
nent is excluded in the DIS case: at large Nf, there 
is no solution for the corresponding Sudakov distribu- 
tion function, if one defines the Sudakov effective cou- 
pling by eq.|0 with T new (u) — ► Tsdg(u) (moreover this 
effective coupling turns out to have renormalons). In 
the Drell-Yan case, the corresponding Sudakov distri- 
bution function does exist, but implies 0(1/Q) power 
corrections (and renormalons in the associated effective 
coupling). With the option to incorporate all constant 
terms into the Sudakov exponent disfavored, there re- 
mains the alternative to absorb into the exponent a 'min- 
imum' set of constant terms, realized |l5j by the (univer- 
sal) choice T new (u) = 1 corresponding to the 'simple' 
coupling eq. (|10|l . A variant of this 'minimal' option, in- 
spired by the work of ^(|, would allow to correlate the 



AC F (H^(Q 2 ) - 2H new (Q 2 )) 



In 



~7WT 



choice of 'left-over' constant terms functions H n 



and 



H^y n °t included in the exponents S new and S^fy to the 
absolute ratio | "^^^ | of time-like and space-like quark 
form factors: one could require (at least at large Nf) that 



Thus, given a 'natural' choice of G^fy , hence of H^y , 
the corresponding DIS functions H new and G new would 
be fixed. An important issue is that of universality 
[l7| (here extended between space-like and time-like pro- 
cesses): given an ansatz for a universal Sudakov effective 
coupling at large Nf (e.g. ea . l(Tfj|) ) . one should check 
(work in progress) whether it remains universal at the 
perturbative level at finite Nf, using the same Sudakov 
distribution functions (e.g. eq.JHJ) and (|f lfl ) found at 
Nf = oo. Further insight could be afforded by a better 
understanding of OPE at large N in the DIS case. What- 
ever the correct choice of the Sudakov distribution func- 
tion turns out to be, it remains for future phenomeno- 
logical work to determine the corresponding form of the 
(possibly non-perturbative) Sudakov coupling at finite 
Nf for each process, and test for eventual deviations from 
universality. This step requires a parametrization of the 
IR part of the coupling, which can be viewed as an alter- 
native to the shape function approach 0, but does not 
involve any explicit IR cut-off. 

I thank M. Beneke, Yu.L. Dokshitzer, J-P. Lansberg, G. 
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cussions. I am indebted to S. Friot for the result quoted 
in [15|. I also wish to thank the referee for thought- 
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